home *** CD-ROM | disk | FTP | other *** search
/ Best of Shareware / Best of PC Windows Shareware 1.0 - Wayzata Technology (7111) (1993).iso / mac / DOS / PROGRAMG / GRAD / DERS2.FOR < prev    next >
Text File  |  1993-01-06  |  1KB  |  33 lines
  1. C ===  Derivating with respect to: 
  2. C     X Y Z                                                                
  3. C
  4.       SUBROUTINE FP (X, Y, Z, F)
  5. C == Fletcher and Powell, 1963 ==
  6.       IMPLICIT REAL*8 (A-H,O-Z)
  7.       PI = 3.14159
  8.       PI2 = 2 * PI
  9.       R_1=(X+X)/2./DSQRT(X*X+Y*Y)
  10.       R_2=(Y+Y)/2./DSQRT(X*X+Y*Y)
  11.       R = DSQRT(X*X + Y*Y)
  12.       IF (X .LT. 0) GO TO 10
  13.       THETA_1=(-Y/X)/X/(1.+(Y/X)**2)/PI2
  14. C*** WARNING: New identifier THETA_1 too long ***
  15.       THETA_2=1./X/(1.+(Y/X)**2)/PI2
  16. C*** WARNING: New identifier THETA_2 too long ***
  17.         THETA = DATAN(Y / X) / PI2
  18.         GO TO 20
  19.    10 CONTINUE
  20.       THETA_1=(-Y/X)/X/(1.+(Y/X)**2)/PI2
  21. C*** WARNING: New identifier THETA_1 too long ***
  22.       THETA_2=1./X/(1.+(Y/X)**2)/PI2
  23. C*** WARNING: New identifier THETA_2 too long ***
  24.         THETA = (PI + DATAN(Y / X)) / PI2
  25.    20 CONTINUE
  26.       F_1=100.*(2*(Z-10*THETA)**(2-1)*(-10*THETA_1)+2*(R-1)**(2-1)*R_1)
  27.       F_2=100.*(2*(Z-10*THETA)**(2-1)*(-10*THETA_2)+2*(R-1)**(2-1)*R_2)
  28.       F_3=100.*2*(Z-10*THETA)**(2-1)+Z+Z
  29.       F = 100. * ((Z - 10 * THETA)**2 + 
  30.      :    (R - 1)**2) + Z*Z
  31.       RETURN
  32.       END
  33.